Effect of the dynamical phases on the nonlinear amplitudes' evolution 
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In this Letter we show how the nonlinear evolution of a resonant triad depends on the special com- 
bination of the modes' phases chosen according to the resonance conditions. This phase combination 
is called dynamical phase. Its evolution is studied for two integrable cases: a triad and a cluster 
formed by two connected triads, using a numerical method which is fully validated by monitoring 
the conserved quantities known analytically. We show that dynamical phases, usually regarded as 
equal to zero or constants, play a substantial role in the dynamics of the clusters. Indeed, some 
effects are (i) to diminish the period of energy exchange r within a cluster by 20% and more; (ii) 
to diminish, at time scale r, the variability of wave energies by 25% and more; (iii) to generate 
a new time scale, T >> t, in which we observe considerable energy exchange within a cluster, as 
well as a periodic behaviour (with period T) in the variability of modes' energies. These findings 
can be applied, for example, to the control of energy input, exchange and output in Tokamaks; for 
explanation of some experimental results; to guide and improve the performance of experiments; to 
interpret the results of numerical simulations, etc. 

PACS numbers: 47.27.Ak, 47.27.ed, 52.25.Fi 



1. Introduction. Nonlinear resonances are ubiqui- 
tous in physics. Euler equations, regarded with various 
boundary conditions and specific values of some param- 
eters, describe an enormous number of nonlinear disper- 
sive wave systems (capillary waves, surface water waves, 
atmospheric planetary waves, drift waves in plasma, etc.) 
all possessing nonlinear resonances Nonlinear reso- 
nances appear in mechanics Q, astronomy Q, medicine 
[if, etc., etc. In this Letter we will regard the simplest 
nonlinear resonant systems corresponding to the 3-wave 
resonance conditions. Examples of these nonlinear reso- 
nant systems are, in order of simplicity, triads (which are 
integrable), and small groups of connected triads which 
are known to be important for various physical applica- 
tions (large-scale motions in the Earth's atmosphere 5|, 
laboratory experiments with gravity-capillary waves 6j, 
etc.). Dynamical system for a triad will be regarded in 
the standard Manley-Rowe form: 

B\ = ZB^B^, £>2 = ZB\B^, B3 = —ZB\Bi, (1) 

where (B\, B2, B3) are complex amplitudes of 3 reso- 
nantly interacting modes Bj expi(kj ■ x — ui(kj)t), while 
the corresponding resonance conditions are 

W (k 1 )+cj(k 2 )-^(k 3 ) = 0, k 1 +k 2 -k 3 = 0, (2) 

where uj(k) is the dispersion relation and k is the 
wavevector. 

Sys. {!]) has been studied both in its real and com- 
plex form by numerous researchers (e.g., 0], [H, etc.). 
If regarded in the amplitude-phase representation Bj = 
\Bj\ expiOj, Sys.([T]) is equivalent to a system for the 3 
real amplitudes \Bj\ and the phase combination ip = 
— the individual phases 9j being slave variables 
and can be obtained by quadratures The dynamical 
equation for ip is also known ([3|, p. 43, Eq.(28)). Still, a 



sort of general misunderstanding persists, concerning the 
relevance of ip for the general dynamics of the system. It 
is a common belief that for an exact resonance to occur, 
it is necessary that the phase tp is either zero ([l(J, p. 132, 
Eq.(6.7); [nj], p.156, Eq.(3. 26.19), etc.) or constant (e.g. 
0). 

The main goal in this Letter is to show that in the case 
of generic initial conditions, the phase ip which we call 
from here on dynamical phase affects the evolution of 
the amplitudes and therefore, has a direct impact on the 
behaviour of any physical system governed by a triad 
as well as small clusters of resonant triads. As it was 
shown in [13] , some of these clusters are described by in- 
tegrable systems and for them a complete set of conser- 
vation laws (CLs) was given explicitly, showing that dy- 
namical phases are relevant in the determination of these 
CLs. In this paper we present differential equations for 
the two independent dynamical phases appearing in the 
butterfly, a resonance cluster formed by two triads con- 
nected via one mode. We investigate the integrable case 
of butterfly by solving numerically the reduced evolution- 
ary differential equations shown in [12| for the phases and 
amplitudes. The numerical integration is fully validated 
by monitoring the conservation laws, known analytically 
from our previous work [l2| ■ We study the effects of the 
phases on the modes' amplitudes, and some physical im- 
plications are briefly discussed. 

2. Triad with complex amplitudes. The CLs for the 

Sys. CO) have the form 



7 13 = |i3 1 | 2 + |i?3| 2 , 
I T = lm(BiB 2 B^) , 



I23 = |-E>2 



(3) 



and their knowledge is enough for the explicit solution 
to be constructed. Here conserved quantities /13 and I23 
are not energy and enstrophy anymore but their linear 
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FIG. 1: Color online. Plots of the modes' amplitudes and dynamical phase as functions of time, for a triad with Z = 1. For each 
frame, <p(t) is (red) solid, Ci(t) is (purple) dotted, C2(t) is (blue) dash-dotted, C?,(t) is (green) dashed. Upper panel. Initial 
condition a — 0.7 and conserved quantities I23 ~ 1, I13 = 1.1 for all frames. Initial conditions for the phase and corresponding 
cubic conserved quantity: Left: tp = 0, It = 0. Middle: ip = 0.1, It = 0.041. Right: tp = tt/2, It = 0.408. Lower panel: 
Conserved quantity /13 = 1.1 for all frames. Initial conditions and corresponding conserved quantities: Left: a = 1.56, ip = 0.01 
and hi = 1,I T = 3.2 x 10~ 5 . Middle: a = 0.7, 93 = 0.01 and I 23 = 0.1, It = 5.1 x 10~ 4 . Right: a = 0.7, ip = tt/2 and 
I23 = 0.1, It = 5.1 x 10 -2 . Here, horizontal axe denotes non-dimensional time; vertical left and right axes denote amplitude 
and phase correspondingly. 



combinations l5j . The analytical solution of Sys.{T]) can 
be found in [8(], as well as the equations on two of the 
three phases 9j (the standard amplitude-phase represen- 
tation Bj = CjCxp(i9j) is used). Below we use slightly 
different notations, introduced in [l2j|, because they are 
more convenient for further studies of bigger groups of 
connected triads. The equation for the dynamical phase 
can be easily deduced and reads as 

^ = -/ T (cr 2 + c 2 - 2 -c 3 - 2 ). (4) 

Combining §5§ and it is easy to see that the con- 
straint It = implies that dynamical phase vanishes, 
i.e. ip = 0. 

If we put It = 0, the solution for amplitudes takes a very 
simple and familiar form: 

'd(i) = dn((-i + t ) z if) 
< C 2 (i) =cn((-i + i ) z^TTz,!^) ^ (5) 
C 3 (t)=sn[(-t + t ) zV7^,^)V^ 

where to is defined by initial conditions and can also be 
written out explicitly. 

3. Triad, case It 7^ 0. Now we present some results for 
the cases when It ^ but otherwise the resonance con- 
ditions in the standard form ^ are satisfied. In FigfTJ 
the characteristic evolution of amplitudes is shown, de- 
pending on the value of It- To characterize the initial 
conditions, the variable a = arctan(C3/C2) has been cho- 
sen for a triad and variables a a = arctan(C3 a /C2 a ) and 
Q'6 = arctan(C3b/C2h) have been chosen for a butterfly. 
These variables appear naturally from the explicit form 



of the corresponding Hamiltonians, and simplify the form 
of the dynamical systems and conservation laws (see [l2[ 
for more details). 

As it is shown in Figfl] when initial dynamical phase is 
zero it will remain zero at all times, but the amplitudes 
will change sign periodically (upper- left panel). When 
dynamical phase is initially very small but non-zero, the 
amplitudes become purely positive and dynamical phase 
will have abrupt jumps, at those times when the ampli- 
tudes used to change sign (lower-left panel). Physically, 
in terms of squares of amplitudes, the dynamics in both 
cases is quite the same (figure not shown). However, the 
phase's dynamics, with its periodic motion, is revealed 
in the non-zero case. As it is shown in upper panel (left, 
middle and right), non-zero dynamical phase influences 
the evolution of amplitudes, so that as initial tp increases 
from to 7r/2, the range of amplitude variations decreases 
from 1 to 0.1 and the period of the motions decreases 
from 5 to 3. 

In Fig[TJ lower panel, we show that the notion of A-mode 
(active) and P-mode (passive) introduced in (com- 
pare to stability criterion [3]) is useful also in the case of 
non-zero dynamical phase (p. A-modc is the mode with 
the highest frequency, and two other modes are called P- 
modes. In the pictures, C\ is P-mode and C3 is A-mode. 
Lower-middle picture: when initial value of amplitude 
C\ » C3, C2 (4 times on figure), the P-mode C\ keeps 
its energy and the A-mode C3 interacts strongly with the 
remaining P-mode C2. If C2 >> Cs,Ci then the situ- 
ation will be qualitatively the same, with the P-mode 
C 2 keeping the energy. Lower-left figure: on the other 



3 




FIG. 2: Color online. Upper panel. Plots of dynamical phases as functions of time, for a butterfly with Z a — Z b = 1. For 
each frame, <p a (t) is (red) solid and <Pb(t) is (black) dashed. Initial conditions a a = 0.3, a;, = 0.7; and conserved quantities 
lab = 1.1, 723a = ^236 = 0.5 for all frames. Initial conditions for the phases and corresponding cubic conserved quantities: 
Left: (p a = ifib = 0,1b = 0. Middle: tp a = 0, <f>b = tt/2,Ib = 0.13. Right: tp a — (pi — ir/2, Ib = 0.36. Middle panel. Plots 
of modes' amplitudes of a butterfly as functions of time, same initial conditions and parameters as in the corresponding left, 
middle and right frames of upper panel. Connecting mode C\ is (blue) dotted, C2a is (green) long-dash-dotted, Csa is (green) 
long-dashed, C26 is (purple) short-dash-dotted, C3& is (purple) short-dashed. Lower panel. Plots of amplitudes squared, as 
functions of time, same initial conditions and parameters as in the corresponding left, middle and right frames of upper and 
middle panels. Colors as above. To facilitate view, Cfb, C\ h are shifted upwards by the value 1, and C\ is shifted upwards by 
the value 1.75. Here, in all three panels, the horizontal axe denotes non-dimensional time; the vertical axe denotes phase in 
the upper panel and amplitudes in the middle and lower panels. 



hand, if C3 >> C\, C 2 , then a completely different time 
evolution is observed and all modes interact. 

One more important general feature of the dynami- 
cal phase is shown on the upper- and lower- right panels. 
Indeed, independently of details of the initial values of 
C\ 1 C 2 , C3, the variation range of the amplitudes is mini- 
mized when the initial condition for the dynamical phase 
ip is equal to n/2. This can be used in real physical sys- 
tems in order to control the exchange of energy between 
resonant modes, at no energy cost: the choice of initial 
dynamical phase does not change the energy of the sys- 
tem, which is a sum of squares of amplitudes, obviously 
independent of the dynamical phase (and of any phase, 
for that matter). 



Butterfly with complex amplitudes. As it was 

131 ] . clusters formed by two triads a and b con- 
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nected via one mode can have one of three types accord- 
ingly to the types of connecting mode in each triad: PP-, 
AP- and AA-butterfly. In this Letter, a PP-butterfly is 
taken as a representative example. Dynamical system 



describing evolution of a PP-butterfly has the form 

B\ = Z a B 2a B 3a + Z b B2 b B 3b , 
-82a = Z a B{B 3a , B 2 b = Z b B{B 3b , 
-63a = —Z a B\Bia , B 3b = —Z b BiB 2b ■ 



(6) 



where notation B\ = B\ a = B\ b is chosen for the 
amplitude of the mode common for both triads while 
£?2a, B 3a , B2b, B 3b are other four modes of the butterfly 
cluster. The set of constructed conservation laws reads 
as 



^23a — \B2a\ 



\B 3a \ 



hsb = I-B2&I + \B3b\ 
1|~ -I- |-D3a|" ■+■ l-°3fc| 

Ib = ^{Z a BiB 2a B 3a + ZbBiB 2 bB^ b ) 



I ab = \B 1 \ 2 + \B 3a \ 2 + \Br ]2 



(7) 



Similar to the triad, the use of standard representation 
Bj = Cj exp(i9j) shows that the Sys.((7|) has effectively 
3 degrees of freedom and two dynamical phases are im- 
portant: lfi a = 6ia + 9 2a - #3a, </5fc = &lb + #26 - #35, 

with the requirement 6\ a = 0\b which corresponds to the 
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choice of connecting mode. Accordingly, equations on 
the dynamical phases take form 



^ = -/ B (cr 2 + c 2 - 2 -c 3 - 2 ), 



¥b = -Ib{C\ + C. 



21, 



°3t> 



(8) 



In order to study the effect of non-zero dynamical phases 
for butterfly, we performed numerical simulations for the 
integrable case Z a = Z^. This allows us to compare the 
results with the triad, which is just a particular case of 
the integrable butterfly. 

In FiglU left column: we show phases, amplitudes 
and amplitudes squared (energies) for initial conditions 
<^a = — 0. The dynamics is quite similar to the triad 
with initial condition tp = shown in Fig[T] upper left. In 
Figf2l middle column: phases, amplitudes and energies 
are shown for initial conditions ip a = n/2,(pb = 0, while 
in the right column the same data are presented, for ini- 
tial conditions tp a = ifb — ir/2. We observe from FigsQ] 
and [2] some effects from the dynamical phases (p a (t) and 
(fb(t) of a butterfly (correspondingly, the phase ip(t) of 
a triad): to diminish the period of energy exchange r 
within a cluster by 20% and more; to reduce the vari- 
ability of wave energies by 25% and more; to generate a 
new time scale, T » r, in which there is considerable 
energy exchange within a cluster, as well as a periodic 
behaviour (with period T) in the variability of modes' 
energies. 

All computations have been done using Mathematica 
and we have validated the code by checking the corre- 
sponding conservation laws, particularly those of cubic 
and quartic dependence on the amplitudes (introduced 
in [l2|]). These conservation laws are stably conserved 
during the whole numerical simulation, within a relative 
error of 10 -12 . 

A comment regarding the ergodicity of the integrable 
butterfly. We observe in a parametric plot of cos(<p ) 
vs. cos(ipb) as functions of time (figure not shown), 
that the seemingly periodic motions are indeed process- 
ing with precession speed depending on the initial con- 
ditions. This is a generic feature of integrable systems 
which are not superintegrable. 

5. Conclusions. Effects of non-zero dynamical phases 
should be taken into account in the following situations. 

• To control energy input, exchange and output in lab- 
oratory experiments, e.g. in Tokamaks. Indeed, a possi- 
bility of concentrating energy in a small set of drift waves 
via some instability mechanisms has been conjectured by 
V.I.Petviashvili some 15 years ago As soon as energy 
is concentrated in a resonant cluster, mode amplitudes 
can become dangerously large. In it was shown that 
the appearance of resonances can be completely avoided 
by special choice of the form of the laboratory facilities 
which, of course, is too costly a game with Tokamaks. 
On the other hand, adjustment of dynamical phases can 
diminish amplitudes of resonantly interacting drift waves 



10 times and more for the same technical equipment. 

• To gain more insight into the phenomenon of zonal 
flows in plasmas which are now regarded as the main 
component in all regimes of drift wave turbulence. "The 
progress of plasma physics induced a paradigm shift from 
the previous 'linear, local and deterministic' view of tur- 
bulent transport to the new 'nonlinear, nonlocal (both 
in real and wave number space), statistical' view of tur- 
bulent transport. Physics of the drift wave-zonal flow 
system is a prototypical example of this evolution in un- 
derstanding the turbulence and structure formation in 
plasmas" ( fl7| ) . In [T3|, a modulational instability of 
Rossby /plasma drift waves leads to generation of zonal 
jets through a process in which the wave amplitudes are 
initially well-approximated by a kite, a cluster consist- 
ing of two triads connected via two modes. The study 
of the behaviour of associated dynamical phases could 
lead to a deeper understanding of zonal jet formation. 
Structure formation (for 3-wave resonance processes) is 
presented in (lij , examples of non-local interactions are 
given in [2(| as well as the cases of 'weak' locality (waves 
with wave numbers of order n and n 2 can form a reso- 
nance cluster, n and n 3 - can not); effects of initial energy 
distribution among the modes of a cluster are studied in 
[l3| . In this Letter we identify dynamical phase as an ad- 
ditional important parameter for any theoretical study of 
nonlinear wave systems. We would like to point out that 
the explicit equation for the dynamical phase of a triad 
has been known for more than 40 years (Eq. (3.31), [2l|) 



in plasma physics and probably even earlier in nonlin- 
ear optics. We consider as our material impact in this 
metier the detailed study of dynamical phases' effect on 
nonlinear evolution of a triad and a butterfly. 

• To guide and improve the performance and analy- 
sis of laboratory experiments. We showed that non-zero 
phases can dramatically reduce the variability of the os- 
cillations (Figs. [1] and [2 left columns). It would then 
be possible to tune initial conditions and/or forcing (in 
rotating water tanks, for example) in order that the mea- 
surement of the resonant modes' oscillations be less sub- 
ject to errors. In Q results of laboratory experiments 
with gravity-capillary waves are presented, and corre- 
sponding dynamical system for three connected triads 
is written out and solved explicitly. The authors report 
qualitative agreement of the observation with the solu- 
tions of the dynamical system though the magnitudes of 
observed amplitudes are higher than those theoretically 
predicted. For all calculations the dynamical phase of 
the initially excited triad was set to ir which might be 
the source of this discrepancy. 

• To interpret the results of numerical simulations. 
For instance, in Q a generic model of intra-seasonal os- 
cillations in the Earth's atmosphere has been presented 
which describes the processes with periods of the order 
of 30-90 days. In this model, dynamical phase has not 
been taken into account yet. A new time scale, T » r, 
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which is clearly observable in Fig. [2j corresponds, for the 
resonant triads of atmospheric planetary waves, to peri- 
ods of the order of 2-5 years. This is the time range of 
climate variability. This means that in numerical mod- 
eling of the climate variability the control of dynamical 
phase is a matter of prodigious importance. 
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